Voltage control of magnetocrystalline anisotropy 
in ferromagnetic — semiconductor/piezoelectric hybrid structures 
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We demonstrate dynamic voltage control of the magnetic anisotropy of a (Ga,Mn)As device 
bonded to a piezoelectric transducer. The application of a uniaxial strain leads to a large re- 
orientation of the magnetic easy axis which is detected by measuring longitudinal and transverse 
anisotropic magnetoresistance coefficients. Calculations based on the mean-field kinetic-exchange 
model of (Ga,Mn)As provide microscopic understanding of the measured effect. Electrically induced 
magnetization switching and detection of unconventional crystalline components of the anisotropic 
magnetoresistance are presented, illustrating the generic utility of the piezo voltage control to provide 
new device functionalities and in the research of micromagnetic and magnetotransport phenomena 
in diluted magnetic semiconductors. 

PACS numbers: 75.47.-m, 75.50.Pp, 75.70.Ak 



The control of magnetic properties of ferromagnetic 
materials by electrical means is an important prerequi- 
site for successful implementation of spintronics in infor- 
mation processing technologies, and for major advance- 
ments in sensor and transducer applications. The ma- 
terial class of (III,Mn)V dilute magnetic semiconductors 
(DMSs) plays a major role in this research because it 
combines ferromagnetic and semiconducting properties 
in one physical system and because of the potential com- 
patibility with modern microelectronic technologies. The 
intense research into DMS systems has led to several suc- 
cessful demonstrations of direct gate control of ferromag- 
netism via carrier redistribution |l|, 0] . However, the high 
carrier-doping levels (>^ 1%) necessary to achieve fer- 
romagnetism render the direct gating efficiency in these 
semiconductor materials relatively low. 

While unfavorable for the direct gating, the large con- 
centrations (Fermi energies) of the spin polarized holes 
that mediate ferromagnetic coupling between the Mn lo- 
cal moments produce large magnetic stiffness, resulting in 
a mean-field like magnetization and macroscopic single- 
domain characteristics of these dilute moment ferromag- 
nets. At the same time, magnetocrystalline anisotropics 
derived from spin-orbit coupling effects in the hole va- 
lence bands remain large, even at these high hole con- 
centrations. This leads to the sensitivity of magnetic 
easy-axes orientations to strains as small as 10~^ 
So far the strain effects have been controlled by lattice 
parameter engineering during growth Q or through 
post growth lithography 0, |3,[3|- 

It has been demonstrated that sizeable strains can be 
induced in GaAs structures using a piezoelectric trans- 
ducer Q. Here we utilize this technique to demon- 
strate the dynamic voltage control via strain of the mag- 
netic anisotropy in a (Ga,Mn)As device bonded to a 



piezo-transduccr. We demonstrate that the favorable mi- 
cromagnetic characteristics of DMSs and the relatively 
simple band structure allow for a microscopic descrip- 
tion of these effects on an unprecedented level of ac- 
curacy c omp ared to metal ferromagnet/piezoelectric de- 
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vices [y|, [iUI, [111, [13 ■ Finally we discuss the realization 
of electrically induced magnetization switching and of 
the detection of unconventional crystalline components 
of the anisotropic magnetoresistance (AMR). These are 
two examples illustrating the generic utility of the piezo 
voltage control to provide new device functionalities and 
in the research of micromagnetic and magnetotransport 
phenomena in DMSs. 

The 25 nm thick Gao.94Mno.06 As epilayer was grown by 
low-temperature molecular-beam-epitaxy on GaAs sub- 
strate and buffer layers (see for details). The mate- 
rial is under compressive in-plane strain of ~ 3 x 10"'^ 
14| due to the lattice mismatch with the GaAs. From 



SQUID magnetometry on an unstrained sample the mag- 
netic easy axis is in-plane in a direction determined 
by competition between biaxial [100]/ [010] and uniaxial 
[lIO] anisotropics. The uniaxial term dominates above 
10 K. At 50 K the cubic and uniaxial anisotropy con- 
stants determined from hard axis magnetisation curves 
are Kc = 85 Jm'^ and Ku = 261 Jm'^ ± 20%. 

A (Ga,Mn)As Hall bar, fabricated by optical lithogra- 
phy, and orientated along the [110] direction, was bonded 
to the PZT piezo-transducer using a two-component 
epoxy after thinning the substrate to 150 ± 10 /im by 
chemical etching. The stressor was slightly misaligned so 
that a positive/negative voltage produces a uniaxial ten- 
sile/compressive strain at « —10° to the [110] direction. 

The induced strain was measured by strain gauges, 
aligned along the [110] and [110] directions, mounted on 
a second piece of 150 ± 10 /im thick wafer bonded to the 
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FIG. 1: The longitudinal resistances, Kxx ((a) and (b)) and 
the transverse resistances Kxy ((c) and (d)) as a function of 
magnetic field for angles close to the easy axes (30° at -150 
V and 80° at +150 V). The curves close to the easy axes in 
each case are relatively flat as a function of field, indicating 
small rotation of the angle of the magnetisation. T=50 K 



piezo-stressor. DifFerential thermal contraction of GaAs 
and PZT on cooling to 50 K produces a measured biax- 
ial, in plane, tensile strain at zero bias of 10^'^ and a 
uniaxial strain estimated to be of the order of ~ 10"'^ 
[ist which could not be accurately measured. At 50 K, 
the magnitude of the additional strain for a piczo voltage 
of ±150 V is approximately 2 x 10~^. 

In this study the orientation of the in-plane mag- 
netisation of the (Ga,Mn)As Hall bar was determined 
from the measured longitudinal and transverse AMR. 
To a good approximation (« 10%), these are given by 
l^PxxIPav = Ccos2(j) and Pxy/pav = Csin20, where 
is the angle between the ma gnet isation direction and the 
Hall bar (current) direction [16[. 

Figure 1 shows magnetoresistance measurements at 
50 K for external magnetic field sweeps at constant field 
angle measured from the Hall bar direction. The 
strongly 0-dcpendent low-field magnetoresistance, which 
saturates at higher field is due to AMR, i.e., to mag- 
netisation rotations. We have subtracted an isotropic, 
^-independent magnetoresistance contribution, approxi- 
mated as linear, from the measured longitudinal resis- 
tances. 

When the external field is close to the magnetic easy 
axis, the measured resistances at saturation and rema- 
nencc should be almost the same and a significant mag- 
netoresistance due to rotation of the magnetisation can 
only be present at very low applied fields. For external 
fields away from the easy axis, large magnetoresistances 
corresponding to large rotations of the magnetisation ori- 
entation are present. This enables us to determine the 
easy axis directions within ±5°. 

The effect of the piezo-stressor is clearly apparent in 



Figure 1. At 50K, SQUID measurements show that the 
magnetic easy axis is oriented along the [110]- direc- 
tion for the as-grown (Ga,Mn)As wafer, consistent with 
\Kc\ < \Ku I . The easy axis for the Hall bar bonded to the 
stressor rotates to an angle </) = 65° upon cooling to 50 K 
due to a uniaxial strain induced by anisotropic thermal 
contraction of the piezo stressor [15[. Application of a 
bias of -1-150 V to the stressor causes the easy axis to ro- 
tate further to (j) = 80° while for -150 V it rotates in the 
opposite sense to = 30°. This directly demonstrates 
electric field control of the magnetic anisotropy in our 
(Ga,Mn)As/PZT hybrid system. 

The magnetic anisotropy for our system can be de- 
scribed phenomenologically by an energy functional 
E{M) = -Kc/Asin^ 2(f) + sit? <P + K'^ sin^(0 + <^o), 
where the last term, with 0o ~ 10° i is due to the mis- 
aligned stressor. The observed behaviour is then con- 
sistent with the (Ga,Mn)As being in tensile strain along 
the axis of the stressor on cool down and applied posi- 
tive (negative) voltage increasing (decreasing) this strain. 
Note that the misalignment allows smooth rotation of a 
single easy axis in the experimentally accessible voltage 
range. 

We now calculate the expected magnetic anisotropy 
characteristics of the studied (Ga,Mn)As/PZT system by 
combining the six-band k • p description of the GaAs 
host valence band with the kinetic-exchange model of 
the coupling to the local Mnca d^-moments 0, Q . This 
approach is well-suited to the description of spin-orbit 
coupling phenomena near the top of the valence band 
whose spectral composition and related symmetries are 
dominated by the p-orbitals of the As sublattice, and also 
provides straightforward means of incorporating lattice 
strains [1, H, [a]- 

Due to the presence of unintentional compensating de- 
fects in (Ga,Mn)As films, the concentrations of ferromag- 
netically ordered Mn local moments and holes cannot be 
accurately controlled during growth or determined post 
growth 17| . We therefore consider in our analysis uncom- 
pensated MuGa moment concentrations within an inter- 
val X = 3 — 5% which safely contain the expected value of 
X. The magnetocrystalline anisotropy constants depend 
strongly on the local moment density and the hole com- 
pensation ratio p/Nmu-, where p is the hole density and 
Nun is the concentration of Mn ions. For fixed p and 
-/Vmii, the cubic term Kc-, calculated without adjustable 
parameters, agrees with the measured 50 K value for 
p/Nun — 0.6 — 0.4 for x = 3 — 5% in good agreement 
with the estimated compensation ratio in our as-grown 
material 

The origin of the uniaxial anisotropy term in bare 
Ga,Mn)As wafers is not known, but it can be modelled 
18| by introducing a shear strain Bint along the [110] 
axis. For p/Nyin = 0.6 — 0.4 we obtain the experimen- 
tal T = 50 K value of for compressive shear strain 
Cint = 3 — 2 X 10"'^ within the considered range of x's. 
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FIG. 2: (a) The microscopic E{M) curves for the three piezo 
voltages. cj> is the angle of the magnetisation with respect to 
the Hall bar. (b) The longitudinal AMR from theory calcula- 
tions with a non saturating magnetic field of 20 mT rotated in 
the plane of the film (c) The experimental AMR curves with 
a field of 40 mT rotated in the plane of the film, pav is the 
pxx averaged over 360° in the low field regime. 9 is the angle 
of the magnetic field with respect to the Hall bar. T=50 K 



The calculations reproduce the measured V easy axis 
for a tensile strain of Cstr = 6 — 4 x 10"'*, along the 
stressor axis and the experimental easy axes for ±150 V 
are obtained by increasing/decreasing the estr strain by 
3 — 2 X 10^'*. These changes in strain agree with the 
measured values for ±150 V, and the V strain due to 
differential contraction is of the expected order. The re- 
sulting microscopic E{M) curves for the three voltages 
are shown in Figure 2(a). 

The magnetoresistance calculated microscopically 



from the same band structure model combined with 
Boltzmann transport theory [31 gives AMR at satura- 
tion of the same sign and comparable magnitude to the 
experiment if we assume the above compensation ratios. 
This allows us to microscopically simulate AMR mea- 
surements assuming the single domain behaviour. In 
Figure 2(b) we show the results of simulations and in 
Figure 2(c) experimental data for the situation where a 
magnetic field of magnitude smaller than the saturation 
field is rotated in the plane of the (Ga,Mn)As epilayer. 
Both theory and experiment show that these AMR traces 
are no longer sinusoidal since the magnetisation does not 
track the applied rotating field. Ranges of magnetic field 
angles 9 for which resistance is more slowly varying cor- 
respond to angles close to the easy axis, where the mag- 
netisation vector lags behind 9. On the other hand ro- 
tation around the hard axis is more abrupt, and in this 
region the AMR can develop hysteretic features whose 
widths increase with decreasing magnitude of the rotat- 
ing field. At +150 V the hard axis is close to the Hall 
bar axis resulting in sharper minima than maxima in the 
corresponding experimental and theoretical AMR traces, 
while the trend is clearly opposite for the -150 V bias 
data, consistent with the easy axis directions obtained 

from the field sween mepisnrpmpnts 
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FIG. 3: (a) Low field magnetic hysteresis curve at -f 150 V. 
The field is swept from saturating negative field at 165° to the 
position show by the black arrow. Then (b) the piezo voltage 
is swept inducing a rotation of the angle of the magnetisation, 
indicated by the red arrows. Numbered arrows represent the 
order and direction of the voltage sweeps. T=30 K. 
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Having established a microscopic understanding 
of the control of the magnetic anisotropy in the 
(Ga,Mn)As/PZT hybrid system we now proceed with the 
demonstration of an electrically induced magnetisation 
switching. The bias-dependent hysteresis loops which al- 
low for such a reversal process are shown in Figure 3(a). 
With the piezo voltage at -I- 150V, the initial magnetisa- 
tion state is prepared by sweeping the external magnetic 
field from negative saturating field at 165° to the position 
shown by the black arrow. This causes the magnetisation 
to rotate from 165° to 260°, at B=0 T (i.e. along the easy 
axis at -f 150 V), then to 275° for the small positive field 
of approximately 18 niT (marked by the black arrow). 
Then, with the external magnetic field held constant the 
piezo voltage is swept (Figure 3(b)) and the magnetisa- 
tion rotates from 275° to 25° (i.e. close to the easy axis 
for -150 V) resulting in a change of R^y as shown by the 
red arrows. This sequence switches the magnetisation 
from the 4th to the 1st quadrant, where it remains for 
subsequent voltage sweeps. The magnetisation can be 
switched back again by reversing the sequence, with the 
magnetic field set to the opposite polarity. 

Finally we report on the detection of an unconven- 
tional crystalline component of the AMR allowed by the 
piezo voltage control. The AMR in (Ga,Mn)As is known 
to consist of a non-crystalline component, reflecting the 
symmetry breaking imposed by a preferred current direc- 
tion, and crystalline terms reflecting the underlying crys- 
tal symmetry. The crystalline terms typically represent 
10% of the total AMR in 25nm (Ga,Mn)As layers [H. 
Figure 4 shows the change in the longitudinal lS.pxx/ Pav 
and transverse l^pxy/ Pav components of the AMR for 
piezo voltages of ihl50 V. /^^pxx — Pxx ~ Pav^ 

and pav is 

the average of pxx over 360° in the plane. The distortion 
of the lattice by the piezo transducer leads to modiflca- 
tion of the crystalline components of the AMR, shown in 
the flgurc by subtracting the curves at piezo voltages of 
±150 V. This modiflcation represents « 10% of the total 
AMR and is comparable to the absolute magnitude of 
the crystalline terms. The appearance of a fourth order 
term in the transverse AMR is expected under a uniaxial 
distortion , but this higher order term was not found 
to be significant in mistraincd 25nm (Ga,Mn)As layers 

To conclude, wc have demonstrated the voltage con- 
trol of the magnetic anisotropy and non-volatile switch- 
ing of the magnetisation direction in (Ga,Mn)As induced 
by strain applied with a piezoelectric transducer. Micro- 
scopic theory calculations capture the physics involved. 
These techniques open up new avenues for exploring a 
variety of micromagnctic and magnetotransport phenom- 
ena in DMS systems. 
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FIG. 4: The change in (a) the longitudinal A.pxx/ pav and (b) 
the transverse l^p^y/pav components of the AMR for piezo 
voltages of ±150 V. 
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